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Faddeev-Yakubovski equations are solved numerically for *He tetramer and trimer states using 
realistic helium-helium interaction models. We describe the properties of ground and excited states, 
^1^' and we discuss with a special emphasis the *He-*He3 low energy scattering. 
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^ ■ I. INTRODUCTION 

Helium atom is one of the simplest fcw-bodysystcms. Nowadays, its structure can be described theoretically with 
an accuracy better than the spectroscopy one [l| . In spite of this simplicity, the compounds of Helium atoms display 
a series of unique microscopic and macroscopic physical phenomena. 
^ \ The closed shell electronic structure, as well as the compactness of the Helium atom, makes it the most chemically 
c/2 ■ inert noble gas. Nevertheless, the very weak van der Waals attraction between two distant He atoms is responsible 
^ \ for the fact that at very low temperatures, both bosonic ^He and fermionic '^He liquify. In addition, the extreme 
^ ■ weakness of the He-He interaction is decisive in explaining why liquid Helium is the only known superfluid. 

As much exciting is the physics of the atomic Helium nanoscopic structures (He multimcrs). Recent inspiring 
experiments 0, 1^ have demonstrated the existence of bound diatomic ''He systems (dimers) - the ground state with 
the weakest binding energy of all naturally existing diatomic molecules. Moreover, diffraction experiment projecting 
a molecular beam of small He clusters on nanostructured transmission gratings, allowed a direct measurement of the 
' dimer average bond length value < R >= 52 ± 4 A. Such a large bond, compared to the effective range of the He-He 
, potential, enables one to estimate its binding energy E « ^ = 1.1 + 0.3/ — 0.2 mK, as well as the '*He-'*He 

l/^ ■ scattering length oq = 104 + 8/ - 18 A. 

The extreme weakness of the *He dimer binding energy requires a precise theoretical description of the ab initio 
0^ ■ He-He potential, which results from subtracting the huge energies of separated atoms. Nevertheless, several very 
] accurate theoretical models 0,|HS0Hil3 were recently constructed, predicting Helium properties in full agreement 
with experiments. These effective He-He potentials are dominated by the strong repulsion (hard-core) at distances 
Rne-He ^ 2 A, where the two He atoms are "overlapping". At larger He-He distances, the weak van-der-Waals 
Q i attraction takes over, creating a shallow attractive pocket with a maximal depth of Vb « 10.9 K, centered at Rne-He ~ 
2>A. The strong repulsion of the He-He potential at short distance, allows to treat accurately atomic He systems without 
having to take into account the internal structure of single He atoms. 
^ '• The physics of small He clusters is an outstanding laboratory for testing different quantum mechanical phenomena. 
The bound state of ^He dimer, at practically zero energy, suggests the possibility of observing an Efimov-like state in 
the triatomic compound The existence of this ^He dimer bound state very close to threshold, is also responsible 
for a resonant '*Hc-'*He S-wave scattering length ao ~104 A, which is more than 30 times larger than the typical length 
scale given by the Van der Waals interaction I ~ 3A. The small He clusters low energy scattering observables should 
5— j ' therefore be little sensitive to the details of He-He interaction and should thus exhibit some universal behavior. This 
sharp separation of different scales can be exploited, making He multimcrs a perfect testground for Effective Field 
Theory (EFT) approaches [H IH IH El . 

We present in this work a rigorous theoretical study of the smallest '*He clusters. There exist in the literature a 
large number of theoretical calculations on triatomic ^He (trimer) 0, ^3 ■ However the four-atomic *He system 
(tetramer), being by an order of magnitude more complex in its numerical treatment, remains practically unexplored. 
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Our study tries to fill up this gap. by providing original calculations for tetramer bound and scattering states. In 
addition, some existing ambiguities in triatomic '^He calculations jl^ are discussed. 

To describe ^He tetramer, Faddeev-Yakubovski (FY) equations in configuration space are solved. In some cases, this 
method may be cumbersome and numerically expensive. It constitutes nevertheless a very general and mathematically 
rigorous tool, with the big advantage over many other techniques, that it enables a systematic treatment of bound 
and continuum states. 

The paper is structured as follows: in section^we present the theoretical framework used; in scction lllll we highlight 
our ^He trimer results. Section Hvl deals with four-atomic helium systems and in section O we discuss the possible 
existence of the rotational states in three-atomic and four-atomic^He systems. Section IVII concludes this work with 
the final remarks. 

To describe the interaction between the helium atoms, we have used the potential developed by Aziz and Slaman Q, 
popularly referred to as LM2M2 potential. There exists several equivalent interaction models. However, they all have 
similar structure and quantitatively provide very close results |l7l Il9l | . 

All results presented in this paper are restricted to the bosonic ^Hc isotope; therefore in the following, we omit the 
mass number 4 and refer to ^He as He. All calculations use i- = 12.12 K-A2 as the input mass of He atoms. 



The Schrodinger equation is the paradigm of non-relativistic quantum mechanics. However this equation is not able 
to separate different rearrangement channels in the asymptote of multiparticle wave function (w.f.). Thus it does not 
provide explicitly a way to implement the physical boundary conditions for scattering w.f., which are necessary to 
obtain a unique solution. Faddeev [2(il | have succeeded to show that these equations can be reformulated by introducing 
some additional physical constraints, what leads to mathematically rigorous and unique solution of the three-body 
scattering problem. Faddeev's pioneering work was followed by Yakubovsky. In 0| , the systematic generalization of 
Faddeev equations for any number of particles was presented. 

One should mention that Schrodinger equation may still be applied when solving few-particle bound state problem. 
However then one must deal directly with total systems w.f., which is fully (anti)symmetric and has quite complicated 
structure. Exploiting the knowledge of systems symmetry one often tries to decompose this w.f. into only partially 
(anti)symmetrized components, which had simpler structure and are more tractable numerically. One practical way 
is to decompose systems w.f. into so called FY components. Then FY equations often have obvious advantage 
over Schrodinger one, since it deals directly only with FY components by avoiding construction of full system wave 
function. 

In that follows, we describe only four-particle FY equations; three-particle Faddeev equations are self-contained in 
four-particle ones. 



II. THEORETICAL FRAMEWORK 



A. Faddeev-Yakubovski equations 
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FIG. 1: Faddeev- Yakubovsky components K and H. Asymptotically, as z ^ oo, components K describe 3-1-1 particle channels, 
whereas components H contain asymptotic states of 2-f 2 channels. 
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The calculations presented in this work were performed using the framework of differential FY equations developed 
by S.P. Merkuriev and S.L. Yakovlev |23,|23- The major step in these equations is the representation of the systems 
w.f. as a sum of 18 components: 12 of type K, describing the asymptotic behavior of different 3+1 particle channels, 
and 6 of type H, describing the systems decomposition into two clusters of two particles: 

* = E^4fc+ E (1) 

i<j i<j,k<l 

Here (ijkl) indicates cyclic permutation of particle indices (1234). These FY components are coupled by 18 FY 
equations. Since all particles are identical, several straightforward symmetry relations can be established between 
different FY components. Using these relations, the number of independent FY equations and components reduces 
to two: 



E-Ho-Vj Kt,^s = V{P+ + P-) [(1 + sP3i)Kt2,3 + Hum] (2) 
E-Ho-v)Hi2M = 1>P [(1 + eP34)i^l\3 + ^12 



12,34 



Here P+ = P12P23, P = P23P12 and P — P13P24 are particle permutation operators. Hq is the kinetic energy 
operator of the system and V is the potential energy operator for the particle pair (12). The coefficient e is a Pauli 
factor: s = 1 for bosonic systems and e = — 1 for systems of identical ferniions. 

By applying combination of permutation operators to the FY components of Eq. (Q, one can express the w.f. of 
the system by means of two non-reducible FY components: 

* = [1 + (1 + P+ + P-)sP3i] (1 + P+ + P-)Kt2,3 + (1 + P+ + P')il + P)Hi2M (3) 

It is convenient to treat Eq. Q using relative reduced coordinates. These coordinates are proportional to well 
known Jacobi coordinates, which are simply scaled by the appropriate mass factors. We use two different sets of 
relative reduced coordinates, defined as follows: 



^7. — > _ /o (™i+™j+"'fc)'"i ri^ _ mj r i+rrij r j+rrik r k \ 

for _ftr— type FY components, mj and ~r*i are respectively the i-th particle mass and position vector. To describe 
i7— type components, another set of coordinates is more appropriate: 



m =J'27^iri-r,,) (5) 



^. . > _ /o (mi+mj)(mk+mi) / mk^k+mi'r'i _ raj r j+nij r j 

Using these coordinates, the kinetic energy operator is expressed as a simple sum of Laplace operators: 

Ho = -t^Ar (A, + Ay + A,) . (6) 

Zm m 

Another big advantage of these coordinate sets is the fact that the degrees of freedom due to the center of mass motion 
are separated. 

The dimension of the problem can be further reduced by using the fact that an isolated system conserves its total 
angular momentum and one of its projections J7z. We deal with systems of bosonic He atoms in their ground state, 
which have a total spin equal to zero. In this case, the system w.f. is independent of the spin and the orbital angular 
momentum is conserved separately. Using this fact, we expand FY components on a partial-wave basis (PWB) of 
orbital angular momentum and, omitting the spin: 



a — (^a; :^ y xy 



(7) 
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In this basis, the kinetic energy operator reads: 



i^o = — 

m 



.2 _ ^^2 _ .2 , lAlx + 1) + 1) ^.^z + 1) 

^ y a;2 y2 ^2 



In Eq. (UI), we have introduced the so caUed FY amphtudes {x,y, z) and {x,y, z), which are continuous 
functions in radial variables x,y and z. The symmetry properties of the w.f. with respect to the exchange of two He 
atoms impose additional constraints. One should have amplitudes only with an even value of the angular momentum 
Ix, whether these amplitudes are derived from K ov H FY components. 

In addition, for _ff — type amplitudes J-^ {x, y, z), the angular momentum ly should be even as well. The total parity 
of the system 11 is given by (_)'^+'!'+'' independently of the coupling scheme {K or H) used. 

By projecting each of the Eq. (O to the PWB of Eq. a system of coupled integro-differential equations is 
obtained. In general, PWB is infinite and one obtains thus an infinite number of coupled equations. This obliges 
us, when solving these equations numerically, to make additional truncations by considering only the most relevant 
amplitudes, namely those which have the smoothest angular dependency (small partial angular momentum values 
Ixyly and Iz). 



B. Boundary conditions 

Equations ((SJ are not complete: they should be supplemented with the appropriate boundary conditions for FY 
components. It is usual to write boundary conditions in the Dirichlet form, which at the origin should mean vanishing 
FY components. However, the existence of a large strong repulsion region (hard-core) corresponding to the inner part 
of He- He potential brings additional complications. The relevant matrix elements from the physical interaction region 
(shallow attractive well) fade away in front of the huge repulsive hard-core terms, thus resulting in severe numerical 
instabilities. 

However, such a strong repulsion at the origin simply indicates that two He atoms cannot get arbitrarily close 
to each other: for a repulsive region of characteristic size r/j, the probability that two particles get closer to each 
other at a certain distance r = c < rh will be vanishingly small. It means that the w.f. of the system vanishes in 
part of a four-particle space, inside six multidimensional surfaces Vij = c, where Vij is the distance between particle 
i and j. The most straightforward way to improve numerical stability would be to avoid calculating the solution in 
at least part of the strong repulsion region, and to impose by hand to the wave function to be equal to zero in this 
region. Nevertheless, due to the complex geometry of this domain in the nine-dimensional space of particle relative 
coordinates, this method is not easy to put in practice. 

A nice way to overcome this difficulty was proposed by Motovilov and Merkuriev [23 | . The authors showed that an 
infinitely repulsive interaction at ^ c, generates boundary conditions for the FY components which can be ensured 
by setting: 

(e~Ho- v) A-f^ 3 = 

; ( for a: < c (8) 



and 

^1^2,3 + (P^ + P') [(1 + £^34)^1^2,3 + ^12.34] = 



Hi2M + P[{^ + eP3A)Kt2-s + Hum] = 



In addition, FY components asymptotic behavior should be conditioned as well. For the bound state problem, the 
w.f. of the system is compact, therefore the regularity conditions can be completed by forcing the amplitudes Ta^^^ 
to vanish at the borders of the hypercube [0, X^ax] x [0, ^max] x [0, ^max], i-e.: 

T^^"\X^-,,^,y,z)=T^'^"\x,Y^-,^^,z)^T^'^"\x,y,Z^,,^) = (10) 

In the case of elastic atom-trimer (1+3) scattering, the asymptotic behavior of the w.f. can be matched by simply 
imposing at the numerical border z = ^max, the solution of the 3N bound state problem for all the quantum numbers, 
corresponding to the open channel a^. It worths reminding that only ii'-type components contribute in describing 
3-1-1 particle channels: 

^aAx,y,Z^i,^) = fa^{x,y) (11) 
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Indeed, below the first inelastic threshold, at large values of z, the solution of Eq. (j^jl factorizes into a He trimer 
ground state w.f. - being solution of 3N Faddeev equations - and a plane wave propagating in z direction with the 
momentum /cq.^ — a/ -p- (Ecm — Ehc^ ) ■ One has: 

Here the functions {x, y) are the Faddeev amplitudes obtained after solving the corresponding He trimer bound 
state problem, whereas Enea is its ground state energy; ni^{ka^z) and ji^ika^z) are regularized Riccati-Bcssel func- 
tions. Equations (O in conjunction with the appropriate boundary conditions define the set of equations to be solved. 
The numerical methods employed will be briefly explained in the next subsection. 

Once the integro-differcntial equations of the scattering problem are solved, one has two different ways to obtain 
the scattering observables. The easier one is to extract the scattering phases directly from the tail of the solution, by 

calculating logarithmic derivative ( — ) of the open channel's K amplitude Ua in the asymptotic region: 

tan (5 = Q^^Kj^ y^^^ ^ ^ (12) 

This result can be independently verified by using an integral representation of the phase shifts 

tan <5 = - ^ J ^i^f J; (fc„„ z) (Vu + 1^24 + V34)-^dV. (13) 

$Q„ is the trimer - composed from the He atoms indexed by 1,2 and 3 - ground state w.f. normalized to unity and 
'I' is normalized according to: 

*(f,, y,;, z,) = $(,123) [j^^ tan(J)n,^ (ka^z)] . (14) 

Detailed discussions on this subject can be found in psi. l2^. 

C. Numerical methods 

In order to solve the set of integro-differential equations - obtained when projecting Eq. and the appropriate 
boundary conditions Eq. H8I11() into a partial wave basis - the components JF" are expanded in terms of piecewise 
Hermite spline basis: 

^"ix,y,z) = ^c'^jkiS]{x)Sk{y)Si{z). 

We use piecewise Hermite polynomials as a spline basis. In this way, the integro-differential equations are converted 
into an equivalent linear algebra problem with unknown spline expansion coefficients to be determined. For bound 
states, the eigenvalue-eigenvector problem reads: 

Ax = EBx, (15) 

where A and B are square matrices, while E and x are respectively the unknown eigenvalue(s) and its eigenvector (s). 
In the case of the elastic scattering problem, a system of linear algebra equations is obtained: 

[A - E,^B] X = h (16) 

where a; is a vector of unknown spline expansion coefficients and h is an inhomogcneous term, generated when 
implementing the boundary conditions Eq. Hll|l . For detailed discussions on the equations and method used to solve 
large scale linear algebra problems, one can refer to [26^ . 

III. TRIMER SCATTERING AND BOUND STATES 

As mentioned in the introduction, trimcr states have been broadly explored in many theoretical works. Hvper spheri- 
cal, variational and Faddecv techniques were used to calculate accurately bound state energies |lalrill9ll27ll28|.l29ll3^ 



6 



as well as to test different He-He interaction models. Nevertheless, we found useful to consider these states as a first 
step, before the more ambitious analysis of He tetramer states could be undertaken. Special emphasis will be at- 
tributed to He-He2 scattering calculations, which are less studied and for which some discrepancies were pointed 
out |l8l Hot . Some arguments will also be developed in favor of considering the first trimer excitation as an Efimov 
state. 

We present in TableQlthe convergence of the He trimer states as a function of the partial- wave basis size. It contains 
results for the ground {B3) and first excited (-B3) state binding energies as well as for the He-He2 scattering length 
(oq^^^^). The basis truncations were made by limiting partial angular momentum values, applying selection criteria 
I'x ^ ly l£ Imax- One Can remark that the convergence is monotonic and quite similar for the Hea bound states and 
for the He-He2 scattering length calculations. These results demonstrate that in order to ensure five digit accuracy, 
the partial wave basis must include FY amplitudes with angular momentum values up to 12. However, for reaching 
a 1% accuracy, Imax — 4 turns out to be enough. 



TABLE I: Convergence of He trimer calculations obtained when increasing partial wave basis. In the three columns are 
respectively presented trimer ground (B3) and excited(B3) state energies in mK, as well as atom-dimer scattering length 
(4^+^)) in A. 



^max 


B3 (mK) 


(mK) 


(A) 





89.01 


2.0093 


155.39 


2 


120.67 


2.2298 


120.95 


4 


125.48 


2.2622 


116.37 


6 


126.20 


2.2669 


115.72 


8 


126.34 


2.2677 


115.61 


10 


126.37 


2.2679 


115.58 


12 


126.39 


2.2680 


115.56 


14 


126.39 


2.2680 


115.56 



In Table are summarized some relevant properties of the trimer ground and excited states. Together with 
their binding energies (B) we give some average quantities like kinetic < T > and potential energy < V > values and 
moments of interparticle distances Xij. Our results are in perfect agreement with other existing calculations [rRl29ll30l| : 
the most accurate of them have been included in this Table for comparison. 



TABLE II: Mean values for He trimer ground and excited states. In this table B, T and V indicates respectively binding, mean 
kinetic and potential energies calculated in mK; Xij stands for interparticle distance. 





Ground state 


Excited state 


B (mK) 
< T > (mK) 
<V > (mK) 


126.39 
1658 
-1785 


126.4 [12 
1660 [17] 
-1787 [17] 


2.268 
122.1 
-124.5 


2.265 [IJ] 
121.9 [17] 
-124.2 [17] 




10.95 


10.96 [17] 


104.3 


102.7 [30] 


< Xij > (A) 

< ^--^ > (A-i) 

< X-/ > (A-2) 


9.612 
0.135 
0.0230 


9.636 [12 
0.0233 [31] 


83.53 
0.0267 
0.00218 


83.08 [IJ] 



The situation is more ambiguous for the He-Ife scattering length ag. Results obtained by Sandhas et al. [ig - 
solving Faddeev equations - and Blume et al. |22| - usin g hy phersphcrical harmonics - agreed with each other on 
a scattering length value Gq = 126 A. However Roudnev [l3|, using also Faddeev equations, found a smaller value 
flo = 115.4 ± 0.1 A. It seems that in Sandhas et al. calculation a relatively small grid was employed, which did 
not allow to disentangle the contribution of virtual dimer break-up component in the asymptotic behavior of the wave 
function. A more recent study pj| | of the same authors provided a value of = 115.5 ± 0.5 A, in full agreement with 
Roudnev's result. 

Our calculations are also very close to this value. When using a numerical grid limited to Umax — 450 A we 
obtain ao = 134 A, still far from the final result. Note that hypcrradial grids employed in [l^ were limited to 
Pmax = \/ + y^=460 A. The results of Table HI have been obtained using a grid with Umax = 950 A, i.e. He-He2 
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FIG. 2: "(Color online)" Convergence of He — [He)2 scattering length obtained when enlarging solution domain of Faddeev 
equations in He — (He)2 separation direction (variable y) with additional discretization points. 




distance rmax = -^Vmax = 823 A, which is large enough to reduce the grid-dependent variations to the fourth 
significant digit. As a function of Vmax, the scattering length varies smoothly (see Figure 12) and converges towards 
the value ao(oo) = 115.2 A, very close to the one given by Roudncv oq = 115.4 ± 0.1 A. 

Whether or not the trinicr first excitation is an Efimov [3| state has been an issue of strong polemics [s^, 113 • In 
fact, when using effective range theory and describing system by zero-range interactions, it is not difficult to show 
the appearance of an Efimov state. However, the problem becomes complicated in calculations with finite range 
potentials. Efimov states accumulate according to a logarithmic law 



N : 



and thus the two body scattering length should be each time increased by a factor ~ e'^ « 23 (or the dimer binding 
energy reduced ^ e^'^ 530 times) to allow the appearance of an additional Efimov state. To handle this property, the 
grids employed in calculations should be extremely large and dense, capable on one hand to accommodate extended 
wave functions, and on the other hand to trace the very weak binding of the third particle to the dimer. Such 
requirements can not be fulfilled in the calculations with realistic interactions. Thus the basic clue in claiming 
excited trimer to be an Efimov state is the fact that this state disappears when the interatomic potential is made less 
attractive [TR 0, ■ Actually, if this potential is multiplied by an enhancement factor 7 > 1 then the following 
effect is observed: first, the difference between the trimer excited state binding energy (i?3) and the dimer one {B2) 
increases with 7. Then, for larger values of 7 the difference {B^ — B2) monotonously decreases and for 7 >1.2 trimer 
excited state moves below the dimer threshold and becomes a virtual one. 

The demonstration of Efimov effect can be accomplished only by showing accumulation of new states when dimer 
binding energy decreases. Such a demonstration has been given in ref. [36j using semi-empirical potential HFD- 
B 0. We would like to remark that the formation of new states can be alternatively demonstrated by studying 
the 7-dependence of the three-body scattering length, without the necessity of solving the bound state problem. In 
scattering calculations, the numerical solution of 3-body equations can be reduced to the interaction domain and the 
wave function extended outside this region using analytical expressions |33 . 

In figure Owe display the behavior of the IIe-He2 scattering lengths, when the He-He potential is multiplied by a 
scaling factor 7 < 1, according to V = jV. In this figure, the Hc-He2 scattering length is plotted as a function of the 
fictive dimer binding energy. One can see that, when decreasing 7 - scattering length decreases. However in absence 
of Efimov states one should expect them increasing, since reducing 7, the dimer target becomes larger. Once 7 is 
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FIG. 3; "(Color online)" The change of the atom-dimer scattering length as a function of the dimer binding energy. 



reduced to w 0.990, the scattering length becomes negative and for values 7 w 0.979 it exhibits a singularity going 
from oq = —00 to +00. This singularity corresponds to the appearance of a new trimer bound state (i.e. second 
excited state in He trimer). This analysis clearly demonstrates that He trimer excited state is an Efimov one. It is 
worth mentioning that for an enhancement factor 7 = 0.979, the He dimer binding energy is only 0.046 mK. 



IV. TETRAMER STATES 



The major aim of this article is to provide a comprehensive analysis of the four-atomic He compound (tetramer). 
The first efforts to describe this system were made already in the late seventies by S. Nakaichi et al. using a 
variational method. Latter on, variational Montecarlo techniques were used by several authors [2^ l39l l40l I^H to 
compute the tetramer ground state. These methods are very powerful in calculating = 0+ bound state properties, 
but are seldom generalized to describe excited states and are not appropriate for scattering problem. 

FY techniques were also used by S. Nakaichi et al. to calculate the tetramer ground state binding energy 
and the He-Hca scattering length. However in order to reduce the - at that time - outmatching numerical costs, 
some important approximations were made. The He-He potential was restricted to S-wave and written as a one-rank 
separable expansion and the same expansion was used to represent the FY amplitudes. These approximations led 
to a tetramer ground state which is underbound by 40% with respect to their own variational result [s^. A recent 
attempt to calculate the He tetramer binding energy using S-wave FY equations was done in |43j . although without 
separable expansion of FY amplitudes. 



TABLE III: Convergence of the He tetramer calculations obtained when increasing the partial wave basis of Eq. (|7|l. The two 

columns represent respectively the tetramer ground state binding energy in mK (B4) the and the atom-trimer scattering length 

max{lx,ly,lz) B4 (mK) '^tf^^'' (A) 

"0 34a8 ^^855 

2 505.9 190.6 

4 548.6 111.6 

6 556.0 105.9 

8 557.7 103.7 
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The FY calculations we present here contain no approximation other than the finite basis set used in the partial 
wave expansion Q. This basis set included amplitudes with internal angular momentum not exceeding a given fixed 
value Imax, i-C. fulfilling the condition max{lx,ly,lz) < Imax- The largest basis we have considered has Imax =8, and 
consist of 180 FY amplitudes, a number by two order of magnitude larger than in preceding calculations. Note, that 
the smallest basis, which is often referred to as S-wave approximation, is obtained by fixing Imax = and requires 
only 2 amplitudes, one of type K and one of type H. The convergence is displayed in Table Hill for the tetramer ground 
state binding energy and He-Hea scattering length. 

The corresponding FY K- amplitudes are displayed in Figure^ as a function of the He-Hea distance r = \f^^- ^^'^ 
can see the different scales involved in the bound state and zero energy scattering wave function. 
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FIG. 4: "(Color online)" Comparison of the functional dependence of K-type FY components in one He atom separation from 
Hea core direction. Single, dashed and dot line curves correspond respectively to tetramer ground, excited state and He-Hea 
zero energy scattering wave functions. 



The convergence of the tetramer calculations is sensibly slower than the one observed for the trimer case (see 
Table P). Such a deterioration is due to the complex structure of the involved FY components. Indeed, each He 
atom pair brings an additional hard-core region; the ensemble of these regions crosses over in the multidimensional 
configuration space and results into a single domain with non trivial geometry. Inside this multidimensional 
domain, the total wave function must vanish by cancelling the contributions of the different FY amplitudes, what can 
be achieved only at the price of increasing its functional complexity. 

The convergence of the He-Hea scattering calculations is also slightly slower than for tetramer ground state. This 
is a consequence of the scattering wave function, which presents a richer structure than the ground state one, as 
can be seen in Figure ^ Tetramer ground state has a rather simple geometry: the four He atoms minimize their 
total energy by forming a tetrahedron. The total energy of the 1+3 scattering state (E w —126 mK) is considerably 
larger than the ground state one {E w —558 mK) and provides more fiexibility to each atom. Furthermore, since 
the scattering length value is much larger than the size of the trimer target, the scattering wave function turns to be 
strongly asymmetric. 

Despite these difficulties, the PW basis we have used enables to reach rather accurate results: the tetramer ground 
state binding energy is converged up to 0.4%, while the final variation of the scattering length does not exceed 2%. 
Our best result for the ground state binding energy of He4, B=558 mK, is in perfect agreement with the B=559(l) 
mK value, provided by several variational Monte-Carlo techniques j29l l39ll4(j| . 



It is interesting to compare the bound state result with the effective field theory (EFT) predictions. It follows 
from this approach that systems governed by large scattering lengths should exhibit universal properties. The wave 
functions of such systems have very large extensions with only a negligible part located inside the interaction region. 
The detailed form of the short-range potential does not matter, since the system probes it only globally. It should 
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be therefore possible to describe its action using only few independent parameters (physical scales). One scale is 
obviously required, to fix the two-body binding energy or the scattering length. Keeping fixed the two-boson binding 
energy, the three bosons system would collapses if the interaction range tends to zero, a collapse known as the Thomas 
effect |4^. This indicates that three-body system is sensitive to an additional scale, which can be determined by fixing 
one three-body observable (for instance the 3-particle binding energy or the 1-1-2 scattering length) )45i l46i | . 

It seems [lj| that the four-boson binding energies remains finite if none of its three-boson subsystems collapse. 
Furthermore four- and three-body binding energies are found to be correlated [T^ . The simplest way to establish such 
a correlation law is by using contact interactions. In this way, the two-body binding energy is fixed by the parameter 
strength of zero-range two-body force, whereas the three-body collapse is avoided and its binding energy is fixed by 
introducing repulsive three-body contact term. Using this model Platter et al. 0| demonstrated that inside quite 
a large domain of dimer-trimer binding energy ratio {B2/B^), the correlation between tetramer and trimer binding 
energies is almost linear. In the same study, a numerical approximation of this correlation law was obtained. 

If we take our dimers {BHe2 = 1.30348 mK) and trimers binding energies as the scales for EFT with contact 
interaction, one obtains a tetramer binding energy of 485 and 491 mK respectively, depending on which trimer energy 
- ground or excited ~ is used. This result is only by 13% smaller than our most accurate result. 

One should remark that the EFT scaling laws are derived using contact interactions, which act only in S-wave. 
Therefore it would be more consistent to compare our values obtained using S-wave approximation {Imax = 0). In 
this case, the EFT formulas give respectively = 330 mK and 325 mK depending on which trimer fixes the scale. 

These values must be compared with our result B^ei — 348.8 mK, i.e. EFT is only by 6% off. 

The good agreement between exact and EFT results, proves the great prediction power of the last approach. EFT 
works surprisingly well even beyond its natural limit of applicability: systems with a size significantly exceeding the 
range of interaction. Note, indeed, that the He tetramer and trimer ground states are rather compact objects, in 
which the interatomic separation is about w 10 A, and therefore comparable with He-He interaction range ^ 3 A. 



TABLE IV: Mean values for He tetramer ground and excited states. In this table B, T and V indicates respectively binding, 
mean kinetic and potential energies calculated in mK; Xij stands for interparticle distance. 





Ground state 


Excited state 


B (mK) 


557.7 


559 [42J 


127.5 


128-130 [14] 


T (mK) 


4107 




1900 




V (mK) 


-4665 


-4850 [42] 


-1913 






8.40 




34.4 




< > (A) 


7.69 


7.76 [40] 


24.8 




< ^r/ > (A-^) 


0.156 




0.088 




<xT/>iA-') 


0.0286 


0.0251 [31] 


0.013 





Let us discuss now the He-Hes scattering results. Nakaichi et al. |42|, using S-wave separable expansion of the 
outdated HFDHE2 potential, obtained a large and negative scattering length ao = —116 A. Using the same HFDHE2 
potential and restricting to lmax=0 amplitudes we have got also a negative, although much larger, scattering length 
value ao w-5600 A. This result is however very unstable, as a consequence of the lmax=0 basis inability to describe 
the non-trivial behavior of the FY components in the hard-core region. If we add max(^y,/z) <4 amplitudes, the 
scattering length reduces to fs-SQS A for the same HFDHE2 interaction model limited to S-wave. Similar calculations 
with the LM2M2 potential give ao ~-450 A. As one can see in Table iHll the full potential must be used in order to 
obtain converged results. The presence of He-He interaction in higher partial waves cardinally changes the physics 
of He tetramer: not only the size but also the sign of the resonant scattering length changes, thus indicating the 
emergence of a new excited state in this compound. 

Another effort to evaluate He-Hca scattering length was made by Blume et al. In their work, effective He — Hcn 
potentials were constructed starting from the same LM2M2 He-He interaction model. Without taking into account 
the particle correlations. Blume et al. provided a positive He-Hca scattering length ao=56.1 A. 

As already mentioned, the large positive scattering length indicates the presence of a tetramer excitation close to 
the trimer ground state threshold. Much physics about this tetramer state can be learned by studying the behavior 
of the zero energy 1-1-3 scattering w.f., or what is even more practical, its FY components. These components are 
only partially symmetrized and have a more transparent asymptotic behavior |23l l47l | . The structure of K-type FY 
component displayed in Figure^ proves that this state is the first tetramer excited state: the corresponding open 
channel FY amplitudes have two nodes in z, the He-Hes separation direction. The first node is situated at ~10 A, i.e. 
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inside the Hea cluster, indicating the presence of a compact ground state. The second node is situated at «103.7 A, 
which coincides with the scattering length value. This coincidence is not accidental: at such He-He3 separations, the 
single He atom is already out of the interaction domain of the Hea cluster. Close to this node, the FY components 
reach the well known linear behavior: 

^aS^^y^z) ^ 'pHe3i^^y)i\l^z - 0.0) (17) 

The factor in front of the z-coordinate is due to mass scaling of Jacobi variables Eq. Q . 

A direct calculation of the He4 excited state represents nowadays a hardly realizable numerical task. This state is 
very weakly bound and its wave function very extended. In order to numerically reproduce it, one is obliged to use a 
very large and dense grid. On the other hand, one should be able to ensure a high accuracy to trace a small binding 
energy difference. For the time being, it constitutes an insurmountable obstacle in the context of the large computing 
power demanding four-body calculations. 

Nevertheless the vicinity of tetramers excited state to the He-Hes continuum makes possible the extraction of its 
binding energy from the scattering results. A bound state is identified with a S-matrix Si{k) pole on the positive 
imaginary axis of the complex momentum plane. Since Si{k) is unitary, its general form close to the pole k = ik^ is: 

S,{k) = p^^e^'^'. (18) 

K — IKq 

The momentum feg is related to the tetramer binding energy measured with respect to trimer ground state threshold 
by ABq = °^ , where iJ,3+i = 3m/4 is the reduced atom-trimer mass. 

On the other hand, the well known effective range expansion can be used to approximate the low energy phase 
shifts: 

A:2'+i.ctg(5 = -- + ir,/c2+o(fc2) (^g) 
ai 2 

Combining relations H18|l and p9(l . one obtains an expression for the bound state momentum fco in terms of the low 
energy parameters. For the He-Hea scattering states with Iz = 0, it reads: 

^rokl - fco + — = (20) 
z ao 




FIG. 5: "(Color online)" Extrapolation of He tetramer excited state binding energy from the He-Hea scattering calculations. 
Low energy phase shift fit is made for k*ctg5 values as a linear function of momentum squared (k^). 
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It follows from Eq. (|20|l that in order to obtain the unknown binding energy it is sufRcient to know the low energy 
coefficients Oq and vq. The scattering length value ag has been already determined and the effective range Tq can be 
extracted by fitting with Eq. 1)19(1 the He-Hea low energies phase shifts. The corresponding extrapolation procedure 
is illustrated in Fig. [3 One can see that inside the considered momentum region, close to fc = | fco | , the effective range 
expansion works perfectly well. The obtained binding energy, ABq = 1.087mK , should not suffer much from higher 
order momentum terms in the expansion 1(19(1 . Nevertheless, the He-Hes effective range is pretty large ro=29.1 A and 
influences significantly the extrapolated binding energy value. By ignoring this term, one would get 



AB, 



(0) 



2^3 



«2 ' 



(21) 



and a binding energy of only ASj^^ = 0.751mK will be obtained, i.e. a value by 30% smaller. 

We finally predict the existence of a = 0+ tetramer excited state with binding energy B'^^^ = ABq + Bne^ = 
127.5 mK. This value compares well with the EFT prediction discussed above, which gives the range B"^^^ g 
[128 — 130] mK, a dispersion due to the fact that the interpolation was done for the total binding energy and not the 
sensibly smaller relative value ABq. 

The validity of the procedure we use to obtain the binding energy of the excited tetramer can be verified in the 
dimer case, for which direct bound state calculations causes no difficulty. The procedure is illustrated in Fig. The 
accurate dimer binding energy is BHe2 = 1-30348 mK, a value very close to ABq = 1.087 mK, and allows to control 
the inaccuracy made when disregarding higher order terms in expansion ((19(1 . By considering only two terms in the 
expansion 1(19(1 we have got BHe2 = 1-3036 mK, only differing in the fifth significant digit from the directly calculated 
value. In addition, the He-He potential effective range extracted by fitting the low energy phase shift is ro=7.337 A. 
This value can be independently calculated applying the Wigner formulae from zero-energy scattering wave function, 
which gives ro=7.326 A, in nice agreement with the extrapolated one. Such an agreement demonstrates the validity 
of our approach. 



CO 
D) 

* 



-0.007 



-0.008 



-0.009 



-0.010- 



k*ctg5 = -1/ag+(ry2)*k'+. 



Ma^ =9.97666*10"' +/- 1 .4616*10"' A"' 
r /2=3.66838 +/- 9.7325*10"' A 



0.0000 







Eg=h'k72n= 1.2063 mK 



k„-1/a„-(r/2)*k\=0 
kg=1. 0371x10"' A"' 
E =1.3036 mK 



E =1.30348 mK 

exact 



0.0004 

\e (A-^) 



0.0008 



FIG. 6: "(Color online)" He dimer binding energy extrapolated from He-He scattering calculations. Low energy phase shift 
k*ctg5 is fitted with a linear function of . Exact k*ctg5(k^) values are indistinguishable from the linear fit. 



One should however mention that the very same procedure is not applicable to the trimcr excited state calculation. 
If done, it would lead to a complex momentum = (1.3 ± 0.8z) • 10^^ A^^. This is a consequence of a dimer binding 
energy almost as small as AB*^^^^ = ^he^ ~ ^He2 ■ In this case, the scattering phase shifts are affected by two open 
thresholds and thus a single channel S-matrix theory is not appropriate. In Hc-Hes scattering, the nearest threshold 
He2-He2 opens only for scattering energies Ec,„ «124 mK and thus is well separated from the energy region of interest 
(- 1 mK). 
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It is interesting to compare the effective ranges for He-He, He-He2 and He-Hcs systems. They are respectively 
ro=7.33, 79.0 and 29.1 A. It is not surprising that the atom-dimer effective range is the largest one: this is a 
consequence of a dimcr ground state which is the most extended of the three considered targets. The atom-trimer 
effective range is more than a third of the atom-dimer one and is significantly larger than the range of the He-He 
potential. This suggest that the trimer ground state has a structure with a sizeable probability to find a single He 
atoms separated by 20-30 A apart from its center. 

We can still use the He-Hes scattering w.f. to obtain a relatively good description of the tetramer excited state. In 
fact, the left hand side of FY equations |21 for the bound and zero-energy He-Hes state differ only by a small energy 
term, which has a little effect on the w.f. inside the interaction region, dominated by large kinetic and potential terms. 
These functions sensibly differ only in the He-Hes asymptotes, where they can be described analytically |23,E3 ^ 
tensor product of a strongly bound trimer in its ground state and a plane wave of the remaining He atom with energy 
E-E/zeg. In practice, the two w.f. differ only by FY amplitudes contributing to the open channels. For zero-energy 
scattering, one has a linearly diverging open channel FY amplitude as described in Eq. (|17() . For bound state these 
amplitudes converge very slowly with a small exponential factor: 

(x, y, z) ~ (bHes (x, y)e-Vlko*z ^22) 



The closed channel FY amplitudes rapidly vanish being shrunken by a relatively large e V " exponent, with 
p = \J + + . 

Using this approximation for the tetramer excited state wave function, we have calculated some of its properties, 
which are summarized in Table ITU 



V. TETRAMER AND TRIMER ROTATIONAL STATES 




FIG. 7: Coupling scheme for three-particle FY component. 

Considering that ^He is a spinless atom, the rotational states of "'He-multimers can be classified according to their 
parity (H) and total orbital angular momentum (L). FY components are very useful to classify these states as well as 
to analyse their properties. As we have mentioned in section^ the Z.^ quantum number must be even (see Fig. CI). It 
follows that for He trimer, the iJ^ = 0~ states are forbidden. Further we classify the He trimer states into those which 
can break into dimcr and a freely propagating He atom, and those in which dimcr states are forbidden by symmetry 
requirements. In the first group, we have states with = (1^,2+, 3^, ...). To be decay stable, these states should 
be more bound than dimer. The other group is formed by the = (1+, 2^, 3+, ...) states, which can only break into 
three helium atoms. 

Clearly, the 1~ is the most promising candidate for a trimer stable rotational state. This state is realized with the 
smallest partial angular momentum, i.e. + ^ 1, and consequently contains the smaller centrifugal terms in the 
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Hamiltonian. The non-existence of this state would imply the non-existence of other states of the same decay group - 
i.e. — (1^, 2+, 3~, ...) - since they involve larger centrifugal energies: Ix + ly ^ 2. The existence of stable trimers 
in the second group is very doubtful; the most favorable would be the 2~ one, for which one has already Ix + ly > 3. 

In a similar way, we can classify the rotational states of He tetramer (see Fig.^. In the first group, we have states 
= (1^,2+, 3^, ..) for which the first decay threshold is the trimer ground state. The most promising state inside 
this group is the 1~ one. For this state, the condition Ix + ly + h ^ ^ must be satisfied. 

In the second group, we have trimer decay stable states = 2~, 3+, ..). These states can neither break into 
trimer-atom pair nor into two dimers: their reference threshold is dimer plus two free atoms. = 0^ state represents 
a special case: it can be broken only into four free atoms. 

The FY formalism we are using allows to calculate multiparticle rotational states with the same ease as zero angular 
momentum ones. However, none of trimer or tetramer rotational states have been found stable. The non-existence 
of weakly bound states can be easily verified using low energy scattering techniques for trimer and tetramer states in 
the first decay groups (two cluster). In addition, the calculated scattering length turns out to be an indicator for the 
strength of the interaction between the scattered clusters and thus tests if ncarthreshold bound or resonant states are 
present. 

For atom-dimer scattering in 1~, one obtains a large positive scattering length (ai)3=114.2 A. If He-He interaction 
is changed, this length scales with the dimer size. If we add a small attractive three-body interaction to force a trimer 
binding without affecting dimer, this scattering length becomes smaller. In fact, when the attractive three-body force 
is weak, this scattering length reduces very slowly. Only when this additional force becomes rather strong, close to 
the critical value binding 1~ trimer, the scattering length falls down to -oo. It crosses a singularity, passing from 
-00 to -|-oo at the value where 1~ trimer bound state appears, and then again stabilizes to large positive value. This 
clearly indicates that no trimer rotational state exist with quantum numbers = (1~, 2"*", 3~...). An even much 
stronger attractive three-body force is required to bind the second group = (1+, 2^,3+, ...) states, thus excluding 
the possible existence of bound or even ncarthreshold resonant He trimer states. 

The case of rotational tetramers is identical to the trimers one. A positive scattering length (01)3 sal2.28 A is 
also obtained for He-Hea scattering with = 1^. This value is sensibly smaller than the one for 1^ atom-dimer 
scattering, which simply results from the scaling with the target size. As in trimers case, one has to apply a strong 
additional attractive force to reduce this scattering length to -00, i.e. to force the 1~ tetramers binding. Tetramers in 
the second decay group, as well as — 0^, seem to be even less favorable: they require an even stronger additional 
force to be bound. Wc conclude therefore that no bound or even ncarthreshold resonant states should exist for 
tetramers with ^ 0+. 

VI. SUMMARY 

In this paper, we have outlined Faddeev-Yakubovski equation formalism in configuration space. It enables a 
consistent description of bound and scattering states in multiparticle systems. This formalism was applied to study the 
lightest {N = 2, 3, 4) systems of He atoms using realistic He-He interaction. We have presented accurate calculations 
for bound He trimer and tetramer states, as well as for low energy atom-dimer and atom-trimer scattering. 

Our main results concern the He tetramer states. We have obtained a tetramer ground state binding energy of 
B = 558 mK, in perfect agreement with the most accurate variational results. 

The first realistic calculation of He-Hca scattering length has been achieved, with the prediction ao=104 A. Such 
a large value indicates the existence of a tetramer excited state close to the trimer ground state threshold. 

Its binding energy can hardly be determined in a direct bound state calculation, due to the difficulties in accom- 
modating a very extended wave function. We have shown that this energy can be obtained by applying an effective 
range expansion to low energy 1-1-3 scattering states. We predict the existence of a = 0"*" He tetramer excited 
state with a binding energy of 127.5 mK, situated 1.09 mK below the trimer ground state. 

Finally, we have studied the possible existence of rotational states in three and four atomic He compounds. It has 
been shown that neither the He trimer nor the tetramer have bound rotational (L^ 7^ 0^) states. The existence of 
corresponding ncarthreshold resonances is also doubtful. 
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